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Abstract 



We use a simple QCD-based model to study particle production in S+Au 
collisions at 200 GeV/n. A requisite consistency is met for the hadronic ob- 
servables (vr*^ and vr^ spectra) while pursuing estimates for e^e~ production. 
Since radiative decays of initially produced hadrons has accounted for only a 
portion of the observed dileptons at the CERN SPS, we search for additional 
mechanisms. By including contributions from prompt secondary hadronic 
scatterings, irp — > vre^e", and adding to tt+vt^ annihilation and hadronic 
decays, the "excess" dilepton signal can possibly be interpreted. 
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I. INTRODUCTION 



Measuring and analyzing electromagnetic radiation from heavy-ion collisions represents a 
significant experimental challenge compared to hadronic signals owing to the relatively small 
cross sections. The additional information they provide certainly justifies the undertaking. 
Hadrons produced in the initial stages of the collision interact on average several times before 
leaving the reaction zone. Consequently, any information embedded in hadronic dynamics 
is completely masked by multiple scatterings. Dileptons are not disturbed by the hadronic 
environment even though they are produced at all stages of the collisions as they have long 
mean free paths. They are dubbed "clean probes" of the collision dynamics. 

Recent results from CERN have brought about a surge of activity in search of quanti- 
tative interpretation. The proton-induced reactions (p+Be and p+Au at 450 GeV) are con- 
sistent with predictions from primary particle production and subsequent radiative and/or 
Dalitz decays suggesting that the e~^e~ yields are fairly well understood. Yet, the heavy-ion 
data (S+Au at 200 GeV/n) show a significant excess as compared to the same model for me- 
son production and electromagnetic decays. When integrated over pair invariant mass up to 
1.5 GeV, the number of electron pairs exceeded the "cocktail" prediction by a factor of 5±2. 
It is clear that two-pion annihilation contributes in the heavy-ion reactions as fireball-like 
features emerge and support copious pion production. Vector dominance arguments would 
naturally lead to extra production around the rho mass. Yet, the excess is most pronounced 
between the two-pion threshold and the rho mass. 

The nature of the enhancement suggests several possibilities. Medium modifications 
resulting in a shifted rho mass could be responsible 0. Along these lines, effects arising 
from a modified pion dispersion relation have been investigated 0. Kapusta, Kharzeev, and 
McLerran suggest that the enhancement may be due to a shift in the t]' mass. If medium 
properties are favorable, the rj' can be considered a Goldstone boson with a small mass which 
implies that the cross section for producing an rj' is larger. This in turn could boost the 
number of dileptons in the region of observed excess . 
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Secondary scattering of pions and other resonances has also been studied focusing on 
the role of the ai through np ^ ai —>■ tt e~^e~ . The contribution was shown to be relevant but 
not sufficient for interpreting the data. We extend the secondary scattering investigation 
in the present calculation by including non-resonance dilepton-producing vrp ^ vr e+e~ 
reactions . We shall organize our paper in the following way. In Sec. |I| we discuss the event 
generator providing the basis of the reaction description and we compare to pion spectra from 



experiment. Then in Sec. |T| we introduce an algorithm for estimating secondary scattering. 
The prompt, or non-resonance, vrp — tt e^e^ reactions are modeled with an effective field 
theory. A brief description of gauge-invariantly introducing strong-interaction form factors 
also appears. In Sec. |^ we discuss normalization and acceptance effects, followed by the 
results and conclusions in Sec. 0. 



II. PRIMARY SCATTERING 

Future collider energies, several thousand GeV per nucleon in the center of mass, probe 
distances much smaller than the nucleon. Models must of course incorporate QCD to de- 
scribe the subnucleonic features. Quarks and gluons then comprise the appropriate degrees 
of freedom for a QCD transport theory. They are propagated through spacetime approx- 
imating the dynamics of collisions to be explored at RHIC and LHC [[f|,^p!0|. Evolution 
continues until soft processes dominate and hadronization occurs. 

Whether or not a quark-gluon plasma can be experimentally detected depends largely on 
the characteristics of the collision in its absence, something we shall call background. In order 
to better quantify this background, simulations without "built-in" plasma formation that 
still assume a QCD description of nucleon scattering must be explored. HIJING, developed 
by Wang and Gyulassy is precisely this type of model and has been used to look at 
multiple minijet production, shadowing and jet quenching in pA and AA collisions. 

The simulation we develop is similar to HIJING. It is based on a simple prescription that 
uses QCD to characterize the individual nucleon-nucleon collisions and uses Glauber-type 
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geometry to determine the scaling. The kinematics of the nucleon-nucleon colhsions are 



handled by PYTHIA and JETSET high energy event-generators using QCD matrix 
elements as well as the Lund fragmentation scheme. A somewhat detailed description of the 
model is outlined below. 

• Initialization of nucleons inside nuclei: The nucleons are positioned randomly inside 
each nucleus according to the size of the nucleus and are given random Fermi momen- 
tum in the x-y plane and are given 2;-momentum proportional to the lab energy (or 
center of mass energy). 

• Number of collisions: The number of collisions is determined geometrically [l^. For 
a proton-nucleus collision, 

n{h) = (Jnn j dx dy dz p[\/W- + z'^) (2.1) 
For nucleus-nucleus collisions, 



NifS) = aNN I dxdydzi dz2PA{^x'^ + Z/^ + ^i^) Pb{.\J x"^ + iy - bf + 2^2^) (2.2) 



Picking scattering partners: Two nucleons are chosen at random from each nuclei 
and are allowed to scatter when, and if, they meet several criteria: First, the two 
nucleons cannot have scattered previously. Second, the nucleons must be within one 
cross-sectional radius of one another in the transverse beam direction. 



xt- XpY + {yt-ypY < JonnItt: (2.3) 



Thirdly, the pair must be moving toward one another in the transverse plane. And 
lastly, the center of mass energy, y^, must be above 6 GeV. This limit is chosen because 
it is on the order of the energy where perturbative QCD is no longer applicable. If the 
two nucleons meet these criteria, they are allowed to scatter. 
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Scattering and Rescattering: PYTHIA chooses partons to participate in the hard scat- 
tering from each nucleon. The partons that are chosen, as well as the momentum 



fraction they carry, are based on known parton distributions . After the individual 
partons have had a hard scattering and are color-connected with the diquarks from 
the remaining nucleon, strings are formed. The kinematics of the fragments from the 
string are determined by JETSET and are stored in a temporary array. Any partonic 
radiation that is not color-connected to either string goes directly into the nucleus- 
nucleus final state. This string is then put back into the nuclei and allowed to rescatter 
as a "wounded" nucleon. The wounded nucleon has the string's momentum while its 
position is updated to halfway between the original nucleons' positions through, 

fe.y0(-.!/2)^(^.^). (2.4) 

• Final State: After all nucleons have experienced their geometrically determined num- 
ber of collisions or they have center of mass energies below the cutoff, particles present 
in each nucleon's temporary array constitute the final state. 

Hadronic observables from the model have been compared against data for several sys- 
tems. Although the model is based on very simple premises, it reproduces the main features 
characterizing hadronic final states. Most notable for the present work, we have made com- 
parisons to S+Au hadronic data as well as dilepton CERN SPS data (see section for 
dilepton data). Our model matches the 7r° spectrum from WA80 [1^, as well as the 7r~s 



from NA35 [|T6|. The two comparisons are shown in Figs. 1 and 2, respectively. The pt 
distribution of neutral pions is slightly above the data at low pr and slightly below the data 
for high Pt pions. This might account for part of the difference in kinematic acceptance 
discussed later in Sec |V[ 
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III. SECONDARY SCATTERING 



Dileptons from pseudoscalars (vr^, 77, 77') and vectors {u, , cf)) produced in the primary 
scattering phase are not enough to account for the S+Au data. Our model also incorporates 
secondary scattering of hadronic resonances. All vrs and ps formed during the primary 
collisions of nucleons will have a chance to scatter amongst themselves before decaying. 
The reactions we consider are of two types, one which produces a resonance that decays to 
dileptons and the other which goes to dileptons directly. 

Of the first type, vr+vr^ ^ p*^ — >■ e"'"e~, vr'^p^ — > 7r^e"'"e~, vr^p^ — > a^^ TT^e~^e~ 
have been included. To accomplish these types of scattering, pions and rhos must of course 
appear in the final state of the model described in the previous section. As the default, 
JETSET automatically decays all hadronic resonances, but it also contains provisions to 
prohibit them. We thus allows neutral pions to scatter from charged rhos when conditions 
are favorable. Technically, the steps involved in secondary scattering are similar to those for 
primary scattering. 

• Number of collisions: The number of collisions is again determined geometrically using 
the appropriate density and cross section. 

• Picking scattering partners: A tt+tt^ or vrp pair is randomly chosen and allowed to 
scatter if (1) the pair has not already scattered, (2) the pair is moving toward one 
another in the transverse plane and (3) the pair is within one cross-sectional radius 
of one another in the transverse beam direction. The cross section for creating a 
resonance is taken to be 



with k being the center-of-mass momentum. The full and partial decay widths for 
p° —>■ TT+TT^ are set to 152 MeV. The full ai decay width is 400 MeV and the partial 



partial 



(3.1) 



fc2(^-m,es)2 + rMiV4 



width for ai 



± 



vrp are energy dependent [Q: 
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ra,^.p = ^|^|k| [2(P. ■ Pp? + m,\m^' + k^)] (3.2) 



Resonance formation and deeay: The kinematics of the resonances are determined from 
the pair of hadrons while JETSET decays the resonance into dileptons using appro- 
priate functions for dT/dM"^ and resulting from analysis of the same Lagrangian 
used to derive Eq. ( |3.2|) . 



Dileptons from secondary scatterings of the resonance type increase the number signifi- 
cantly in the region around the p° mass, but not in the region with the largest gap, 0.2 GeV 
<M^< 0.5 GeV. 

The nonresonant component is estimated here by computing the sole process ir^p^ 
7r^e~^e~ and then assuming the other isospin channels contribute equally. The other channels 
involve Feynman graphs that result in a singularity and must be regulated in a full T-matrix 
or some other effective approach ||18|. To this level of estimate, isospin averaging and ignoring 
interference effects between these and the resonant ai contributions is not worrisome. The 
prescription for directly scattering pions and rhos is very similar to the one used for resonance 
scattering as detailed below. 

• Picking scattering partners: The cross section for vr^'p^ vr^p*^ determines how many 
and how often the charged rhos scatter with pions. Using an effective Lagrangian 
approach to be described momentarily, we generate a set of graphs. The resultant 
cross section can be nicely parameterized for ^/s < 2.5 GeV by 

8 / mb \ 

cT(v^) = 1.8mb + — mbGeV^ + 0.5s - — ^ (3.3) 
s2 \GeVy 

• Scattering: Since this is a non-resonant process, the Monte Carlo directly determines 
the kinematics of the final state. Necessary ingredients for such procedures include 
an interaction Lagrangian and a resulting squared matrix element. The Lagrangian 



employed is |19 
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\D^<^\^ - m.2|$|2 _ ]^p^^p,u ^ ]^„,^^py ^-F^^F^^ (3.4) 



where — dfj_ — ieAn — igpp^ is the covariant derivative, $ is the complex charged pion 
field, p^y is the rho field-strength tensor and F^^, is the photon field strength tensor. 
From this Lagrangian, the matrix elements can be determined. In the calculation, the 
graphs involving the oi are neglected as the contribution from ai has already taken 
into account in the resonance portion of the model. There are three graphs whose 
matrix elements are listed below. 

•^1 = ^^2/'^' ,, e^(p„)(2p,+p„)^(2pi + q),u{p.)^''v{p+) (3.5) 
M, = ^^"z^f^P-^^^P^ - ^)M+Pi),^{P-)Yv{p+) (3.6) 

M, = ^e'^{Pa)[X,MP-)Yv{p+) (3.7) 
where X^^ = ag^^ + b(pi^bu + PhuPiJ) + c(pb^pb^ + Pi^iJ)- 

In the t-channel matrix element, a form factor, h^{t) = (A^ — m^)/{h? — t), appears 
to account for the finite size of the mesons. Its presence breaks gauge invariance. In 
order to completely restore gauge invariance, other terms must be added to the four 
point diagram J\4s- a = —1 b = c = (/i+(t) — 1)/(P6 • ? + Pi • ?)• The parameters A 
and m are set to 1 GeV and nip, respectively. 

The absolute square \A4i + M.2 + and da/dM^ were used to Monte Carlo the 

three-body 7r^e~e"'" final state. 

These lepton pairs are of nonresonant origin and are now added to the pairs from reso- 
nance decays. 

Since the pions and rhos are scattering inside the reaction zone, their dynamics are altered 
by the medium. Being of bremsstrahlung type, these mechanisms are therefore susceptible 



to the Landau- Pomeranchuk-Migdal effect [^]. Pions and rhos involved in secondary scat- 
tering will undergo frequent multiple scatterings, and not only with other pions and rhos. 
Therefore, the number of dileptons produced by this scattering is reduced. The reduction 
factor is dependent at minimum on the invariant mass of the lepton pair as well as the mean 
free path of the pions and rhos. We use a reduction 1 — e"*''^'^, where M is the invariant 
mass of the lepton pair and A is the mean free path of the hadrons. For our purposes and 
level of estimation here, we set A to some average value ^ 1 fm p2|. 



The total dilepton yield from our model is the sum of lepton pairs from primary plus 
secondary scattering. The invariant mass distributions of the dileptons from all contributions 
will be discussed in the last section. 

IV. NORMALIZATION 

In order to keep computation time low, the code was run to look at dileptons from tt*^, 
?7, r/', oj, p°, (j) separately. In each run, the hadron considered was allowed to decay only 
into the dilepton channel, all other modes were prohibited. Technically, this is merely a way 
to maximize statistics. To reinstate absolute normalization, all lepton pairs counted were 
multiplied by the branching ratio for the process from which they came. This approximation 
is valid because the lepton decay mode is a comparatively rare event. We have successfully 
used this perturbative technique before to calculate high-energy photon production in the 
framework of BUU transport theory p3|. The secondary scattering resonance production 



was handled much the same way. All resonances produced from secondary scattering decayed 
exclusively into their lepton channels and were later multiplied by the appropriate branching 
ratio. Normalization procedures for non-resonance vrp scattering is somewhat different. 
Instead of a branching ratio from the particle data book as before [^, the fraction of 
dilepton events is based on the calculated cross section for vrp vre+e" divided by the total 
cross section for vrp scattering. For simplicity this total cross section is the sum of the elastic 
cross section, (T^Op±^^±pO plus the resonance cross section, cr^Op±^(j^± (incoherently). 
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V. RESULTS 



The invariant mass spectra of dileptons from the primary scattering part of our model 
for three different systems are displayed in Fig. 3. The top two plots display lepton pairs 
from p+Be collisions and p+Au collisions, respectively. The lower plot shows dileptons in 
question from S+Au collisions. The simulation agrees with the proton-induced data and it 
is reassuring that our S+Au model-results are consistent with the cocktail from the CERES 
collaboration ||T|. Plotting against the actual S+Au data reveals a significant enhancement 
over predictions in the invariant mass region between 200 and 500 MeV. There is also a 
modest enhancement for masses above this range. 

Dileptons from secondary scattering for the S+Au system in our model are shown in 
Fig. 4. The contribution from pion annihilation increases the distribution significantly 
in the rho mass region, but still leaves an excess below the rho mass. We should stress 
that we have taken vacuum rho properties throughout. Radiative Oi decay contributes a 
minimal amount in the excess (or deficit) region, but the contribution from non-resonance 
vrp scattering provides the most significant increase. 

With the inclusion of the secondary scattering previously described, the invariant mass 
distributions of dileptons are shown in Fig. 5. The dilepton spectra from the proton- 
induced interactions are not significantly changed. This result is as expected — dileptons 
from the smaller systems are quantitatively described by primary hadronic decays. The 
proton-nucleus collisions do not create a heated nuclear medium large enough or dense 
enough to bring about significant collective effects. Conversely, the S+Au collision has a 
marked increase in lepton-pair production between an invariant mass of 200 and 500 MeV 
as well as a noticeable increase in the higher mass region. It is not surprising that secondary 
scattering becomes important in the S+Au system, as a dense nuclear medium is created 
during the collision. 

Allowing for the short-comings of our model, results still suggest that secondary scat- 
tering is a viable explanation for the excess found in dilepton data. Inclusion of secondary 



10 



scattering, 1) preserves the consistency the primary scattering in our model has with proton- 
induced data and, 2) enhances the number of dileptons within the region of excess discovered 
in S+Au data. Although this agreement cannot rule out other possible explanations for the 
excess electrons, our model's simplicity is attractive. 
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FIGURES 

FIG. 1. 7r° distribution from WA80 as compared with the model. 

FIG. 2. 7r~ rapidity distribution from NA35 compared to the model. 

FIG. 3. Dilepton invariant mass spectra from primary scattering in the model compared to 
CERES data. 

FIG. 4. Contributions from secondary scattering in S+Au collision. 

FIG. 5. Total dilepton invariant mass distributions, including primary and secondary scattering 
in the model as compared with CERES data. 
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